Regular dendritic patterns induced by non-local time-periodic forcing 
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The dynamic response of dendritic solidification to spatially homogeneous time-periodic forcing 
has been studied. Phase-field calculations performed in two dimensions (2D) and experiments on 
thin (quasi 2D) liquid crystal layers show that the frequency of dendritic side-branching can be 
tuned by oscillatory pressure or heating. The sensitivity of this phenomenon to the relevant param- 
eters, the frequency and amplitude of the modulation, the initial undercooling and the anisotropies 
of the interfacial free energy and molecule attachment kinetics, has been explored. It has been 
demonstrated that besides the side-branching mode synchronous with external forcing as emerg- 
ing from the linear Wentzel-Kramers-Brillouin analysis, modes that oscillate with higher harmonic 
frequencies are also present with perceptible amplitudes. 



I. INTRODUCTION 



Complex patterns observed in nature attract a con- 
siderable interest recently The complicated spa- 
tiotemporal behavior that leads to the formation of such 
patterns is usually associated with the instabilities of sys- 
tems transforming under non-equilibrium conditions. A 
spectacular example of practical importance is dendritic 
growth that appears in anisotropic systems where the 
propagation of the transformation front is coupled with 
the diffusion of a conserved quantity. The respective dif- 
fusional instabilities lead to the formation of a travel- 
ing quasi-steady-state tip that emits surface undulations 
which evolve into side arms in directions determined by 
the anisotropy of the system |4|. Many of the techno- 
logically important materials form by dendritic solidifi- 
cation 0. Analo gous phenomena have been reported 
in biological systems |^ , anisotropic Hele-Shaw cells , 
and in cosmology/particle physics Although exper- 
iments on the freezing of transparent liquids ||9|-p^ clari- 
fied many of the essential features of dendrite formation, 
important questions remained open. For example, the 
role played by thermal fluctuations in side-branching is 
the subject of continuing investigations P, p^ , p^ . It is 
anticipated that a selective amplification of thermal noise 
is responsible for the side-branching and for the observed 
irregularity of dendritic patterns po|-p^ 19 1 . 

While the steady-state behavior of dendritic growth is 
understood fairly well in the framework of recent theories 
(microscopic solvability and phase-field theories 



28 1), less is known of the dynamic response of the 



dendritic morphology to time-dependent external pertur- 
bations. The non-linear systems often develop regular 
patterns under periodic forcing Considering the in- 
herent non-linearity of the equations describing dendritic 
growth, it is reasonable to expect that such periodic per- 
turbations lead to resonance patterns that might be used 



to influence the growth morphology. Besides its scien- 
tific interest, a detailed understanding of such phenom- 
ena could open novel routes for designing materials for 
specific applications. 

Time-periodic forcing of dendritic side-branching has 
already been realized via heating locally the dendrite tip 
by laser beam p9-pl|, or by exposing it to an oscilla- 
tory fiow field |32|. Both methods lead to the forma- 
tion of fairly regular morphologies in a suitable frequency 
range. However, these methods cannot be easily used for 
spatially homogeneous controlling of growth morpholo- 
gies in large volumes. To circumvent this problem, we 
demonstrated [ p3[ | the possibility for tuning the dendritic 
morphology by spatially homogeneous time periodic forc- 
ing via modulated pressure and heating, that change 
the undercooling instantaneously and homogeneously in 
large volumes. 

In this paper we present a detailed analysis of the dy- 
namic response of dendritic solidification to oscillating 
pressure and heating. In identifying the resonance condi- 
tions and the most important process parameters, we rely 
on phase-field modeling, one of the most potent meth- 
ods for describing morphology evolution. The theoret- 
ical predictions will be confronted with experiments on 
quasi-2D liquid crystal layers, known as suitable model 
materials ^ , ^5| -|3^ . The rest of this paper is structured 
as follows. In section II we briefly summarize the phase- 
field equations that incorporate external forcing, and in- 
troduce quantities for characterizing growth morphology. 
Section III is devoted to the experimental setup and the 
details of the measurement techniques. In section IV, 
we present the phase-field predictions and compare them 
with the experimental results. In section V, our findings 
are discussed in the light of theory and experiment on 
tip-localized forcing. Finally, a few concluding remarks 
will be made in section VI. 
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II. PHASE-FIELD THEORY WITH EXTERNAL 
FORCING 

We investigate the growth of the crystalhne phase into 
an undercooled single component liquid at a temperature 
between the melting point and the hypercooling limit 
("unit undercooling") p8[ . In this temperature range, 
the thermal diffusion controls the growth process, as op- 
posed to the (molecule attachment) kinetics controlled 
mechanism taking place below the hypercooling temper- 
ature . Owing to the heat release during solidification, 
the crystalline phase heats up to its melting point, while 
the rest of the latent heat is transported into the cold 
liquid. Under such conditions, the planar front is inher- 
ently unstable against thermal perturbations despite the 
stabilizing effect of the interfacial free energy | . 

The phase-field theory is a powerful tool to study 
moving boundary problems, ft represents a dynamic 
extension of the Cahn-Hilliard theory of first-order 
phase transformations, in which the evolution of a non- 
conserved order-parameter is coupled to thermal or 
chemical diffusion 



Here we use a thermody- 
namically consistent version (Model 1 of |^^) which 
ensures locally positive entropy production and was mod- 
ified to incorporate anisotropics of the interfacial free en- 
ergy and the kinetic coefficient The local state of 
the matter is represented by the phase-field (/ii(r, t) 

To incorporate external forcing, the equations [ ^5y41| 
that describe the evolution of the phase-field and the tem- 
perature field T(r, t) in 2D 

r{e)^ = Q(r)y(</)) - G'(</)) ~ A j 1 [«(0)]2|v0p\ (1) 



{c, + Mel,) - ^]L'{T)}^ + L{T)p'{cl>)^ = kV'T, (2) 

have to be modified. Here r is an empirical relaxation 
coefficient, whose inverse is an intrinsic interfacial mobil- 
ity, and stands for variation with respect to (f>. Other 

notations are as follows: p((/)) — 0^(10 — Idcj) + 60^), 
K^/2 the coefficient of the square-gradient term in the 
entropy functional, that depends on the orientation rep- 
resented by the angle 9 (defined as tan{9) = dy(t)/dx4>), 
G{4)) — (/>^(1 — (/))^/4e, e is a constant that determines 
the height of the intermediate maximum of the double 
well-potential, Q{T) = /J^ ^dr/ = AG(r) 
the Gibbs free energy difference between the liquid and 
solid, while Cp, L{T), and k are the specific heat of the liq- 
uid, the enthalpy difference between the liquid and the 
solid, and the thermal conductivity, respectively. The 
primed quantities denote derivatives with respect to the 
arguments. Note that (f) and p are dimensionless. The 
units of k and are W m^^ K^^ and J ni^^ I^^^, re- 
spectively; while Q, G, Cp and 1/e are measured in J m"'^ 
K-i. 



The dominant effect of the pressure modulation p{t) = 
Po + Ap(t) is a variation of the melting point de- 
scribed by the Clausius-Clapeyron law Tm{p) = T',„(po) + 
l^p/SV / Srm where is the volume change upon solidi- 
fication, and Sm is the entropy of fusion. This enters into 
Eq. (f ) via the Gibbs free energy difference which may 
be approximated as AG{T,p) « Sm[Tm{p) — T] at small 
undercoolings {T„i — T ^ Tm)- The other relevant quan- 
tities [specific heat, heat of fusion, and entropy of fusion 
L'(r)] are far less sensitive to the variation of pressure. 
For example, integrating the Clausius-Clapeyron rela- 
tionship one obtains AL — ApAV[l + Acp/Sm+T{aiVi — 



asVs)/AV], where Acp 



Cp.l 



is the variation of 



the specific heat upon melting, ai^s are the thermal ex- 
pansion coefficients, while subscripts / and s refer to the 
liquid and the solid respectively. Under typical condi- 
tions of our experiments [pressure amplitude ~ 2 bar, 
undercooling of f K, and physical properties of CCH3 
(Appendix A)], the undercooling varies by ^ 6%, while 
the relative change of L (and Sm) is only ^ 10^"^. We 
found that such pressure effects on L and Sm do not in- 
fluence the results perceptibly; i.e. pressure modulation 
needs to be taken into account only via Q{T) in Eq. (1). 
Modulated heating, in turn, appears as a time dependent 
source term in Eq. (2). 

At this stage it is advantageous to introduce dimen- 
sionless variables. This is done similarly to that in [^ . 
However, to avoid complications emerging from an oscil- 
lating melting point, the reduced temperature has been 
redefined in terms of a constant reference temperature 
Tr as u{r,t) = (T - Tr)/iTr - Too), where Too is the 
initial temperature of the undercooled liquid. The math- 
ematical problem is invariant to the choice of Tr so far as 
A = Um — {Tm — Tr)/{Tr — T^o) is inserted as in Eq. (3). 
(With Tr = Tm, the usual form of the phase-field model 
is recovered.) 

The relevant physical properties are combined into 
dimensionless parameters A = Cp(Tr — Too)/L, a = 
V^ujSmL/{12cpao), To = SmOPo/cTo, and e = S/lo, 
a — da/dO, where D is the thermal diffusion coefficient, 
d the interface thickness, and to a reference length com- 
parable to the size of a well developed crystallite. The 
anisotropics of the interfacial free energy a and the ki- 
netic coefficient /3 are represented by the dimensionless 
functions a = 1 + (JnCoslnO) /2 and j3 = \ + f3nCos{n9)/2 
of n-fold symmetry, that are related to the dimensional 
quantities via relationships a{9) = aoCf{9) and (}{9) = 
(3of3{9). Here an and Pn are the anisotropy parameters. 
Lengths and time are scaled by uj and lj^/D, respectively. 
The coefficient k and the constant e are related to S and 
the interfacial free energy p^ ]. 

With these notations Eqs. (1) and (2) transform to 
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- - + 30eaA[u - A(i)]0(l - 0)} 
dx dy dy dx 



|^ + l3O<^^(l-0f^ = V2. + m (4) 

where the modulated pressure and heating are incorpo- 
rated via terms A{t) and B{t), respectively. 

In this work Eqs. (|^) - (|^) have been solved numer- 
ically on an X rectangular grid, A^ ~ 1000, that 
corresponds to an area I y, I oi dimensionless linear size 
I = 5 and a grid spacing Aa; — 0.005. An explicit finite 
difference scheme has been employed in the case of Eq. 
(^, while Eq. has been handled by the alternating- 
direction implicit method which is unconditionally stable 
At t = a crystaUine nucleus (0 = 0) is placed at 
the center of the model area filled by uniformly under- 
cooled liquid (0 = 1, u = —1). The reduced temper- 
ature and the phase-field were kept constant (u = — 1 
and (j) — 1) aX the boundaries of the model space. To 
save CPU time, we used a primitive adaptive mesh tech- 
nique: far from the solidification front Eq. (||) has been 
solved on a rough grid (of spacing 10 • Ax). To model the 
effect of thermal fluctuations and to suppress undesir- 
able lattice effects that favor specific growth directions, 
a spatially and temporally uncorrelated noise of ampli- 
tude 0.01 and zero mean value has been added to the 
dimensionless temperature u in every time step. If not 
stated otherwise, we use the dimensionless parameters: 
a = 350, To = 20, e = 0.005, Ax = 0.005, and time step 
At = lO^''. Owing to the known limitations of phase- 
field modeling |4p^, we performed our calculations at 
relatively large undercoolings, A = 0.40 — 0.68. 

To characterize the spatiotemporal behavior of the 
evolving dendritic morphology the following quantities 
were determined: 

1. We measured the width w of the dendrite behind 
the tip at a distance of ^ = 0.75 (150 pixels). To investi- 
gate this feature for longer times, we performed separate 
simulations on a 1200 x 600 grid, oriented so that the 
larger dimension be parallel with the direction of growth 
(here we used At — % ■ 10^^). The amplitude Ai, in 
the Fourier transform W{f) of w{t) characterizes the re- 
sponse synchronous with external forcing. 

2. The symmetry of the growth patterns was charac- 
terized by the quantities a and 6: 



M 



K 



(5) 



(6) 



where M and K are the numbers of pixels in which 
0.4 < (j) < 0.6 for the whole system and for the lower 
left quarter, respectively. Accordingly, a > measures 
the symmetry of the whole domain with respect to the 
3^— axis, while 6 > quantifies the axisymmetry of a main 
branch. Both parameters are zero for symmetric pat- 
terns. 



III. EXPERIMENTAL SETUP AND CONDITIONS 

The nematic-smectic B (A^ — Sb) phase transition of 
liquid crystals is recognized as an appropriate model of 
crystallization in liquids [^-^ . A specialty of this phase 
transformation is a large anisotropy of the interfacial free 
energy [Q, a property expected to be advantageous for 
the regularization of dendritic morphology on theoretical 
grounds (see section IV. A and |Q). 

To test the predictions, experiments have been per- 
formed on thin CCII3 liquid crystal layers. Ready- 
made cells of E.H.C. Co. (Japan) KSRP-10 (of thick- 
ness 10 /im) and KSRP-02 (2 /im) have been filled with 
CCH3 (Merck, Darmstadt). The surface treatment of 
the bounding glass plates assured the planar alignment 
of both A^ and Sb phases [the directors n(A^) and n{SB) 
that describe the average orientation of the elongated 
molecules are in the plane of the cell], and the conduct- 
ing layers on the bounding plates were used as electrodes. 

For pressure modulation the liquid crystal cell was 
placed into a brass box (see Fig. [^) surrounded by a tem- 
perature controlled hot-stage of accuracy ± 3 mK. The 
gas pressure in the brass box has been regulated by a 
computer controlled solenoid valve system that switches 
on and off an excess pressure Pc preset between to 2 
bar with an accuracy of ± 0.03 bar. This allows square 
wave-like pressure modulations in the frequency range up 
to i> 2 Hz. 

The modulated heat release in the bulk has been real- 
ized by periodically transmitting a high frequency (600 
kHz) electric current through the liquid crystal layer pro- 
duced by gating the signal of a function generator with 
V ^ 1 Hz. The local off-plane heat transport (a precon- 
dition of regulation with oscillatory heating, see section 

IV. A and |Q) is ensured by the quasi-2D sample geome- 
try and by the heat transport through the bounding glass 
plates. 

The growth patterns were monitored in transmitting 
mode via a polarizing microscope equipped with a CCD 
camera, a method relying on the different optical prop- 
erties of the nematic and smcctic B phases. The images 
were stored and processed by a PC. The spatial and time 
resolutions of the system were 512 x 512 pixels and 0.04 s, 
respectively. The calibration procedure with 6.3 x objec- 
tive and 3.2 x projector combination gave scale factors 
of 1.35 ± 0.01 /im/pixel in x direction and 0.95 ± 0.01 
/im/pixel in y direction. 
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The same setup has been used to measure the pressure 
coefficient of the temperature of transformation between 
the smectic-B and nematic phases (see Appendix A). 

IV. RESULTS 
A. Theoretical predictions 

To find resonance patterns we investigated the param- 
eter space defined by the frequency and amphtude of 
the modulation, the undercoohng and the anisotropics 
of the interfacial free energy and the kinetic coefficient. 
The calculations were performed with four- and six-fold 
anisotropics (n — A and 6 in expressions for a and /3). 

Before presenting our results on non-local periodic 
forcing we characterize the unmodulated reference state. 

1. Free dendritic growth 

Independently of the symmetry of the system, a weak 
uncorrelated (noise- induced) side-branching occurs in the 
unmodulated reference states shown in Figs, ^(a) and 
|(a). The width w{t) of a main dendrite arm measured 
at a dimensionless distance of C = 0.75 behind the tip 
[Fig. |4|(a)] indicates surface undulations that can be de- 
composed into a spectrum of oscillations which covers 
the dimensionless frequency range of ~ 30 to ~ 170, cen- 
tered around a characteristic frequency of /o ~ 100 [see 
insert in Fig. ||(b)]. This behavior accords with that 
seen in the phase-field simulations of Karma and Rappel 
[Q , and with experimental results [|l6|,0 . In agreement 
with [|l^, we find that far from the tip the characteristic 
wave length Aq — w//o depends much more weakly on the 
distance from the tip than predicted by the microscopic 
solvability theory |2^,0 , a behavior interpreted in terms 
of the "stretching" of the perturbations [|l8| that travel 
along the curved perimeter of the dendrite. 

Considering that the applied uncorrelated noise probes 
the dynamic response of the system to a broad variety of 
frequencies, it is reasonable to seek "resonance" in the 
vicinity of the characteristic frequency /o of spontaneous 
side branching. 

2. Pressure modulations 

According to the Clausius-Clapeyron law, the modu- 
lated pressure translates into a time dependent melting 
point, and is represented by inserting A{t) = Um{t) = 
aog{t) and B{t) — into Eqs. (3) and (4). Since our 
experimental setup allows only square- wave modulations 
(characterized by the filling coefficient ^ = ton /to, where 
to is the period of oscillations and ton the pulse length). 



the simulations have been performed for this wave form 
(unless stated otherwise). 

(a) The effect of forcing frequency: The frequency de- 
pendence of the growth patterns is shown in Figs, ^(b)- 
(f) and ||(b)-(f). At low forcing frequencies {v < 15), 
the system alternates between two steady-state growth 
modes yielding uncorrelated side-branching that corre- 
spond to the actual undercooling. It is, however, re- 
markable that the switching transient between the con- 
stant pressure stages initiates more pronounced side arms 
than formed during the constant pressure periods. Reg- 
ular growth morphologies [Figs. ||(b)-(e) and ^(b)-(e)] 
are observed in the lower half of dimensionless frequency 
range of spontaneous thermal side branching shown in 
the insert of Fig. |4|(b)]. At frequencies larger than these, 
the formation of side branches cannot follow the exter- 
nal forcing, and the uncorrelated thermal side-branching 
reappears [Figs. ||(f) and ||(f)]. 

In agreement with our earlier results [|3^ , when regular 
morphologies are formed, the side-branching and the tip 
velocity correlate with the pressure modulation. 

The parameters w(t), a, and b reflect the formation of 
regular morphologies. For example, a periodic variation 
of the width w{t) of the dendrite is seen [Fig. ^(a)], which 
correlates with the external sinusoidal forcing and leads 
to the formation of regular side branches. Note that be- 
sides the forcing frequency (that appears with a far larger 
amplitude in W{f) than the spontaneous undulations), 
its second and third harmonics {2i', Sv) are also present 
albeit with an amplitude that diminishes for higher order 
harmonics [Fig. ^(b)]. 

Apparently, periodic forcing with fixed frequency ex- 
cites several surface modulations; one synchronous with 
the forcing frequency (v) and others oscillating with the 
higher harmonic frequencies {2i/,3i>, ...), which indicates 
a dynamic coupling among these modes associated with 
the non-linearity of the governing equations. 

The relative amplitudes of higher harmonics vary with 
the wave form of the pressure modulation yielding dif- 
ferent growth patterns (see Fig. ^). Fairly similar pat- 
terns [cf. Figs. ^ (a) and (b)] are observed for sinusoidal 
and square-wave modulations, provided that for the lat- 
ter ^ = 0.5 [i.e., frequencies {2n + l)^, n = 1,2,... are 
present in the forcing spectrum]. When the filling coef- 
ficient ^ deviates from 0.5, the frequency /phase content 
of forcing varies (e.g., the second harmonic 2iy appears) 
leading to such effects as the formation of side-branches 
of twice the number corresponding to the base frequency 
[cf. Fig. H(a)-(c)]. To understand this behavior we an- 
alyze the frequency dependence of the dynamic response 
of the system to "asymmetric" forcing (^ = 0.3). 

The Fourier amplitudes A^, and A21, (corresponding to 
modes that oscillate with the forcing frequency and its 
double, respectively) and the symmetry parameters a and 
b are presented in Fig. ^ as a function of the forcing fre- 
quency v. Ai, and A21, show maxima if v or 2i/ are close to 
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fexc ~ 70 [Fig. 1(a)], where fexc is the frequency where 
the synchronous response has the maximum amphtude. 
Note that fexc is somewhat below the characteristic fre- 
quency /o ~ 100 of the unperturbed dendritic growth. A 
possible explanation of this deviation might be a different 
"stretching" of perturbations, related to differences seen 
between shapes of dendrite tips formed in free growth 
and under external forcing. It is also remarkable, that at 
low forcing frequencies the amplitude A21, of the second 
harmonic becomes larger than A^. This inversion of the 
relative magnitudes of A^, and A21, is responsible for the 
side-branch "doubling" shown in Fig. |^. A similar phe- 
nomenon is seen in Fig. ^(b), where the second harmonic 
is the dominant mode as indicated by the striking sim- 
ilarity of patterns corresponding to = 20 and 40 [cf. 
Figs. H(b) and (c)]. Note the lack of alternating vari- 
ation of the length of side-branches that appears when 
the synchronous mode is also present with a significant 
amplitude [as in Fig. ||(c)]. 

We find that the amplitude A^ increases approximately 
exponentially with the distance C from the tip of the den- 
drite. In contrast, A21, shows a maximum as a function 
ofC (Fig. 0). 

Remarkably, the symmetry parameters a and b display 
minima at lower forcing frequencies [around v ^ 20; see 
Fig. |^(b)] than the maximum of the Fourier-amplitude 
Ai, . This difference in the positions of the extrema is a 
manifestation of the fact that axial symmetry and fre- 
quency content arc different constituents of regularity. It 
appears that A^, gives a closer representation of the reg- 
ularity recognized by the human eye than the symmetry 
parameters. 

The overall frequency dependence of the pattern for- 
mation supports our anticipation that regularization via 
periodic forcing with fixed frequencies is possible when 
the latter fall in the characteristic frequency range of 
noise-induced side branching. 

(b) The effect of the amplitude of forcing: The Fourier 
amplitudes A,y and A21, and the symmetry parameters 
show that with increasing forcing amplitude ag the reg- 
ularity of the growth patterns increases (Fig. ||). This is 
manifested in monotonously increasing A,y and A21, that 
saturate for large oq. At the same time, a and b de- 
crease with gq , reflecting the increasing symmetry of the 
patterns. 

(c) The effect of undercooling: The response to a given 
pressure modulation correlates with the tendency for 
"natural" side-branching (Fig. ^. While at large under- 
coolings (e.g., A = 0.65), the formation of side-branches 
is rather intense, leading to well developed side-arms 
(both with or without forcing), at A = 0.5 or below free 
growth produces essentially no side-branches, and exter- 
nal forcing is needed to trigger them, also reflected in 
their smaller size. Note the difference in spatial periodic- 
ity, that originates from the temperature dependence of 
the tip velocity. 



The variation of the Fourier amplitudes A^, and A21, 
with undercooling (A) is shown in Fig. |l^. The maxi- 
mum of An might be understood as follows. We found 
that with increasing undercooling fexc increases. Accord- 
ingly, a fixed forcing frequency ly produces the maximum 
amplitude at the undercooling for which fexc ^ i^- For 
the same reason, the maximum amplitude for the second 
harmonic 2^ appears at a larger undercooling. 

(d) The effect of anisotropy: Theoretical considera- 
tions and numerical simulations show that anisotropy 
plays a central role in the formation of dendritic mor- 
phology. It is, therefore, reasonable to expect that it 
has a similar importance in the formation of regular pat- 
terns. This expectation is verified by the dependence of 
growth forms on the magnitude of anisotropies in the 
system (Fig. |ll|). The growth patterns and variations 
of the Fourier amplitudes A^, and and of the sym- 
metry parameters (Fig. ^2|) indicate that the regulariz- 
ing effect of external forcing decreases dramatically with 
vanishing anisotropy, i.e. without well defined orienta- 
tional preferences the external perturbations are unable 
to drive the amoebae-like growth forms into a regular 
pattern. Note that the tip velocity decreases with de- 
creasing anisotropies (note the difference of simulation 
times the snapshots correspond to), the excitation with 
1^ — 60 generates periodic side-arms with increasing spac- 
ing. 

In these calculations we assumed that both the interfa- 
cial free energy and the kinetic coefficient are anisotropic. 
For the sake of simplicity, their ratio was usually kept 
constant Pi/a^ = —2. Our further investigations show, 
however, that resonance patterns form if any of these 
anisotropies (that of a or (3) is sufficiently large. 

3. Modulated heating 

(a) Alternating heating and cooling: They lead to es- 
sentially the same type of resonance patterns as pressure 
oscillations (Fig. |l3|), provided that that the net heat 
production in a period is negligible. A remarkable dif- 
ference is, however, that while the average tip velocity 
vq is essentially independent of the amplitude of pressure 
modulations, it decreases with the amplitude of heat pro- 
duction/extraction (Fig. [1^ ); a phenomenon that might 
be associated with different efficiencies of heating and 
cooling. 

(b) Oscillatory heating: The introduction of a local 
off-plane thermal transport described by B{t) — bog{t) + 
h[u{r,t) — Uoo] was necessary to prevent the melting of 
the crystal for the heating amplitudes needed to gener- 
ate regular patterns. Here ft. < is a dimensionless heat 
transfer coefficient that imitates heat transfer perpendic- 
ular to the plane of the liquid crystal layer towards the 
surrounding of reduced temperature Uoo ■ Since this term 
serves as a local heat sink, the in-plane thermal diffusion 
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becomes less restrictive. This results in the formation of 
more compact objects (see Fig. [l5| ) with thick main arms 
and less developed side-branches. While with an increas- 
ing amplitude of the heat pulses the side-arms become 
more regular, they become smaller due to the dissipated 
heat. 

Much like pressure modulations, oscillatory heating 
leads to an oscillating tip velocity. Owing to a decreasing 
average undercooling accompanied with increasing heat- 
ing amplitude 5o, the period averaged velocity vq of the 
tip decreases roughly linearly with bo (Fig. |l^). At the 
same time, the amplitude of the oscillatory part of the 
velocity Vi, increases about linearly. 

In summary, our investigations imply that under well 
defined conditions both types of non-local forcing can be 
used to control dendritic growth. 

B. Experimental results 



formed at similar reduced frequencies {v/ fo) are com- 
pared [cf. Figs. |(c) and 0(b) of v/fo = 0.40 and 0.42; 
or Figs. |(d) and 0(d) of £>//o = 0.6 and 0.65]. An 
interesting observation is that far from the dendrite tip 
irregularities develop in the lengths of the side-branches, 
despite the even distance of their trunks. We observed 
an analogous phenomenon driven by the interaction with 
neighboring side-arms via thermal diffusion fields in long- 
time phase-field simulations | |45[ |. 

The effect of the forcing amplitude is shown in Fig. 
Pressure oscillations of amplitude of 1 bar effi- 
ciently regularize the side-branch formation. With in- 
creasing amplitude the phenomenon becomes more pro- 
nounced as predicted by the phase-field calculations. At 
the largest pressure amplitudes even the secondary side- 
branches correlate with the pressure oscillations [see Fig. 
^(c)]. Note again that the competition of neighboring 
side-branches leads to irregularities in their lengths. 



In defining the experimental conditions, we utilized the 
results of the computer simulations. For CCH3, the time 
and spatial resolution of our experimental setup is opti- 
mal at an undercooling of 1 K. Comparable variations 
of the undercooling (~ 0.1 ~ 0.2 K) are achievable by 
applying modulation amplitudes as high as 2 bar or 3 
xlO""* W/cm^, as estimated on the basis of the pres- 
sure coefficient of the equilibrium temperature for the 
nematic-smectic B phases (see Appendix A) and other 
relevant properties. 

The appropriate frequency range of forcing has been 
determined by Fourier analysis of the width w{t) of a 
freely growing dendrite measured at 66.5 /im behind the 
tip. Without perturbation, the side-branching is essen- 
tially random [Figs. 0(a) and[l^(a)]. The characteristic 
frequency, identified as the peak of the broad spectrum in 
the insert of Fig. 0, is roughly 1.8 Hz. Accordingly, our 
experimental setup was designed to cover the frequency 
range of - 2 Hz. 



1. Pressure modulations 

The experimental results for pattern formation in the 
2 /im thick liquid crystal cell are summarized in Fig. 0. 
Resonance patterns of fairly regular side-branches are 
observed [Figs. 
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(b)-(f)] at forcing frequencies com- 
parable with the characteristic frequency of the unper- 
turbed state V ~ 1.8 Hz. The figure shows germs that 
nucleated heterogeneously at the same site with essen- 
tially the same orientation several times, yielding repro- 
ducible patterns. Although the respective dimensionless 
undercooling (A = 0.06) is considerably smaller than in 
the simulations (A = 0.4 — 0.68), the observed behav- 
ior follows closely the predictions of the phase-field the- 
ory. This similarity is especially striking when patterns 



2. Modulated heating 

In analogy to pressure modulations, the periodic heat- 
ing experiments on the 10 /im thick liquid crystal cell 
also reveal regular side-branching [see Figs. |l8|(a)-(c)]. 
The formation of side-branches correlates with the exter- 
nal forcing as illustrated in Fig. |l8|(b), where the black 
lines denote the position of the tip at the centers of the 
heating pulses. The correlation is also evident from the 
power spectrum of the width of the dendrite measured 
66.5 ^m behind the tip (Fig. |l|). In fuU accord with 
the phase-field simulations for asymmetric square waves 
[Fig. ^(b)], the power spectrum indicates the presence of 
modes that are either synchronous with the forcing fre- 
quency v = 0.46 Hz, or oscillate with doubled frequency 
(2i>). Even a peak corresponding to the third harmonic 
(3i>) may be identified, although with an amplitude that 
is close to the experimental uncertainty. 

Remarkably, in the response to "symmetric" forcing 
^ = 0.5, the 2v mode is also present (see the Fourier- 
spectrum and the short side-arms of doubled frequency 
in Fig. 0), although this frequency is absent from the 
forcing spectrum. This finding confirms the non-linear 
behavior revealed by phase-field simulations, that higher 
harmonics missing from the forcing spectrum are also ex- 
cited [see W{f) for sinusoidal forcing in Fig. ^. 

Increasing the heating amplitude or ^ so that the 
period-averaged heating power reaches P ~ 10^^ 
W/cm^, the formation of the side-branches is suppressed 
[see Fig. 0(d)], a phenomenon resembling that seen in 
the numerical simulations (Fig. 0). Further increase of 
the heating amplitude (and power) melts the dendrites 
back. 

In line with our theoretical predictions (see Fig. 2. 
of [01), a weak oscillation of the tip velocity has been 
observed that correlates with the forcing. However, its 



6 



amplitude is just above the resolution of the present ex- 
perimental setup. 

A remarkable feature of the "regularized" dendrites, 
not seen in the simulations, is a shift in the position of 
the side-branches on the two sides of the main tip [see 
e.g. Figs. |l^(b) and (c)]. As a result, one cannot use 
the symmetry parameters a and b [defined by eqs. (3) 
and (4)] for characterization of the pattern's regularity. 
For the same reason, the Fourier spectra shown in Figs. 
[l9| and have been evaluated from the half-width of 
the dendrite. Despite the shift of the side-branches, the 
power spectra on the two sides are fairly similar. 

The shift in the position of the side-branches might 
be attributed to the asymmetry of the dendrite tip (due 
to the angular dependence (7{9) [35]), shown by a high 
resolution snapshot of the tip region of a dendrite grown 
freely at AT = 1.0 °C (Fig. ||). Another remarkable 
feature is that the left side of the tip is faceted. The 
first observable surface undulation (marked with arrow 
in Fig. p^ ) that evolves later into a side-arm appears on 
the opposite side. In contrast, the faceted side remains 
smooth up to the same distance behind the tip. 

Finally, one should mention that the electric heating in 
liquid crystals may have side-effects that are not incorpo- 
rated into our phase-field model. Switching the electric 
field on, the orientation of the nematic director n(A^) 
changes from planar to homeotropic (perpendicular to 
the bounding plates), an effect that influences the mag- 
nitude and anisotropy of the interfacial free energy, and 
may induce local flow in the sample. However, these side- 
effects are of minor importance, since they are present in 
full strength much below the electric fleld needed for reg- 
ularization. 



V. DISCUSSION 

In this section we confront our results with those on 
the dynamic response of dendritic growth to local forc- 
ing with flxed frequency as emerging from the micro- 
scopic solvability theory (MST) and experiment. It is 
appropriate to mention, that some features of the MST 
derivation prevent a quantitative comparison. For ex- 
ample, the MST calculations were performed for a 2D 
symmetric dendrite, while the anisotropy, the stability 
coefficient, and the Peclet-number were assumed to be 
small; conditions that are not met in our simulations 
and experiments. Therefore, a qualitative comparison 
is only meaningful. An important further difference be- 
tween previous work and ours is that we used non-local 
forcing, as opposed with the tip-localized forcing assumed 
in the MST. 

Let us first recall some of the MST predictions on noise 
amplification (based on the Wentzel-Kramers-Brillouin 
approximation) that have been tested by comparison 



with numerical simulations [ ^J46| and experiments ]3C| ] 
for fixed frequency perturbations localized at the tip. 

(i) Periodic forcing by a fixed frequency v leads to sur- 
face oscillations of amplitude that increases exponentially 
with the distance from the dendrite tip up to a critical 
distance C,c proportional with \/v'^. For larger distances 
the amplitude decreases and eventually dies away. If, in 
turn, the distance from the tip is fixed and v is varied, a 
peak is observed in the amplitude. 

(ii) The localized wave packets behave differently; they 
grow exponentially as they move to arbitrarily large dis- 
tances from the tip, while the respective characteristic 
wave length increases. 

Our simulations are in line with the previous MST and 
numerical predictions for tip-localized forcing in the fol- 
lowing respects: 

(a) In the vicinity of the dendrite tip, the amplitude 
of the mode synchronous with forcing increases roughly 
exponentially with the distance C from the tip. Owing 
to the excessive computation time and memory needed, 
we were unable to study the decay of this mode in de- 
tail. Nevertheless, as expected from the MST, at large 
frequencies [Fig. ^(f) and ||(f)] the synchronous mode 
can only be recognized in the vicinity of the dendrite tip. 
Another sign, that accords with the presence of a crit- 
ical distance Cc beyond which the excited mode decays, 
is a maximum observed in the amplitude A21, of the sec- 
ond harmonic (Fig. ^ in simulations where the forcing 
spectrum contains with a signiflcant amplitude. 

(b) The amplitude A^, measured at a fixed distance be- 
hind the dendrite tip shows a maximum as a function of 
the forcing frequency [Fig. |^(a)]. 

An interesting feature, revealed by our phase-field 
simulations, is that besides the synchronous mode, the 
Fourier spectrum of the dendrite width contains the sec- 
ond and third harmonics with perceptible amplitudes. 
This finding is confirmed by our experiments. We be- 
lieve that this is the first direct demonstration of such 
non-linear effects in connection with diffusional instabil- 
ity induced dendrite formation during a first-order phase 
transformation. It appears, that this non-linear behav- 
ior is not a peculiarity of spatially homogeneous forcing. 
Although not mentioned in the original works, traces 
of higher harmonics seem to be present in the respec- 
tive power spectra for local forcing as well [^,0. It 
is worth mentioning, furthermore, that a similar behav- 
ior has been observed in the case of anomalous Saffman- 
Taylor fingering [^^, suggesting that the appearance of 
higher harmonic modes in response to periodic exter- 
nal forcing is a common feature of dendrite formation 
whether driven by diffusional or mechanical instabilities. 

Summarizing, in the present state of affairs it seems 
that the only specific feature of pattern formation un- 
der non-local forcing is a simultaneous triggering of side- 
branches (of first, second, and higher orders) on indepen- 
dent crystallites throughout the sample. 
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VI. SUMMARY 

Our computer simulations and experiments demon- 
strated that the dendritic morphology can be regular- 
ized by non-local time-periodic forcing realized by mod- 
ulated pressure and Joule heating. These conditions 
lead to an oscillatory velocity of the dendrite tip, and 
yield side-branches at regular distances, provided that 
the frequency of the modulation is close to the "natu- 
ral" frequency of free side-branch formation. The dy- 
namic response of the system to such non-local forcing 
can be understood in general on the basis of previous the- 
oretical and experimental results on tip-localized forcing. 
However, the non-linear effects, such as the appearance 
of higher harmonic modes, warrant further theoretical 
work. 



(molar mass), p = 895.4 kg/m^ (density at 80.6 °C 
I!), A7J = 6247 J/mol, /Wn^sb/V = -2.6 %, and 
a Is! — 25x 10^"' 1/K (volumetric expansion coefficient in 
the nematic phase |^^) one obtains = 0.033 K/bar 
in excellent agreement with the experimental value. 

Although this value is much larger than for metals, it 
is not unusual for liquid crystals. For example, ~ 0.03 
K/bar has been reported for the nematic - crystal tran- 
sition in PAA, and ~ 0.1 K/bar for the nematic smectic 
transition in p-methoxybenzoic acid [^0|| . Comparable, 
or even larger coefficients have been measured on other 
organic substances such as camphene 0.214 K/atm, pi- 
valic acid 0.0674 K/atm, and succinonitrile 0.0245 K/atm 
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APPENDIX A: PRESSURE DEPENDENCE OF 
THE PHASE TRANSITION TEMPERATURE 

The pressure dependence of the N — Sb phase tran- 
sition temperature has been measured by a procedure 
similar to that described in For different Pe a sin- 
gle Sb monodomain has been kept at constant size by 
controlling the temperature for several hours (until the 
thermodynamic equilibrium state of the system has ap- 
proached). Fig. |2^ shows the change of the phase tran- 
sition temperature depending on the pressure. Linear fit 
on the data gives a slope ^ = (0.032 ± 0.003) K/bar 
for the Clapeyron coefficient. 

The Clapeyron coefficient can be also calculated from 
the relation: 



dr„ tI/Wn-,sb 



dp 



AH 



(Al) 



where T^, AVn^sb and AH are the phase transition 
temperature at atmospheric pressure, the molar volume 
change on transition and the molar latent heat of fu- 
sion, respectively. From the relevant material parame- 
ters of CCH3, namely = 329.45 K pi, m = 233 
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FIG. 1. The experimental setup. 
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FIG. 2. Dendritic patterns of four-fold symmetry pre- 
dicted by the phase-field theory for A = 0.6, a a = —0.12, 
Pa = 0.24 and t = 0.22: (a) Under constant pressure 
yielding A = Um = 0.048; and (b)-(f) in the presence of 
square-wave pressure modulations of amplitude ao = 0.16 
and filling coefficient £, = 0.3. The gray level in the liq- 
uid phase indicates the temperature field. The white line 
marks the position, where the width of the dendrite arm 
was measured. 



FIG. 3. Dendritic patterns of six-fold symmetry pre- 
dicted by the phase-field theory for A = 0.6, ae = —0.04, 
f3e = 0.08 and t = 0.22: (a) Under constant pressure yield- 
ing A = 0.048; and (b)-(f) in the presence of square-wave 
pressure modulations of amplitude ao = 0.16 and filling 
coefficient ^ = 0.3. 
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FIG. 4. (a) The width of the dendrite at a dis- 

tance of = 0.75 behind the dendrite tip without periodic 
perturbation [heavy sohd hne; Fig. |^(a)], in the presence 
of sinusoidal (light solid line) and square-wave [dotted line; 
Fig. 0(d)]. (b) The respective Fourier spectra W{f). 




FIG. 5. Comparison of growth patterns correspond- 
ing to (a) sinusoidal, (b) symmetric (^ — 0.5) and (c) 
asymmetric (^ — 0.3) square-wave forcing of equal period 
lengths {v — 30), and period averaged amplitudes. Note 
the side-branch "doubhng" in panel (c). (Other relevant 
parameters are as for Fig. H.) 
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FIG. 6. The amplitude of the Fourier peaks A^, and 
A21, measured at = 0.75 behind the tip, and the sym- 
metry parameters a and 6 as a function of the driving 
frequency v. Other relevant parameters are A — 0.55, 
^ = 0.3, 0-4 = -0.12 and /34 = 0.24, while ao = 0.14 (a), 
and ao ~ 0.2 (b). 
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FIG. 7. The Fourier amplitudes Ai, (circles) and A21, 
(triangles) as a function of the distance ( from the dendrite 
tip at two undercoolings A — 0.55 (full symbols) and 0.6 
(empty symbols). Note the roughly exponential behavior 
of Ai, and the saturation/maximum in A2u- Other relevant 
parameters are 1/ = 60, ^ — 0.3, ao = 0.16, 0-4 = —0.12 and 
(34 = 0.24. 
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FIG. 8. The amplitude of the Fourier peaks A^, and 
A2,^, and the symmetry parameters a and 6 as a function 
of the forcing amplitude ao for the system shown in Fig. 
id). 
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FIG. 10. The amplitude of the Fourier peaks A^, A2v 
as a function of undercooling A. Other relevant parameters 
are i/ = 60, ^ = 0.3, ao = 0.16, a a = -0.12 and (Sa = 0.24. 




FIG. 11. Pattern formation under periodic forcing as 
a function of anisotropy: 0-4 = —0.08 (a), 04, = —0.04 
(b), 0-4 = -0.02 (c) and = -0.01 (d). Other relevant 
parameters are A — 0.55, v — &Q, = 0.3, ao = 0.16 
and Pa — —2(74. The average growth rate varies with 
anisotropy, therefore, different times are needed to reach 
comparable sizes: t = 0.24 (a), t = 0.3 (b), t = 0.38 (c), 
and t = 0.44 (d). 




FIG. 9. The effect of undercooling on pattern forma- 
tion under periodic forcing, (a) A = 0.65; (b) A = 0.5. 
Other parameters as for Fig. ^(d) except that t = 0.2 and 
t = 0.283 for panels (a) and (b), respectively. 



0.12 
0.10 

0.08 

< 0.06 

< 0.04 
0.02 



0.00 
0.00 



•A, 




a 






1- b 


- 


• 




- A 


• • • 










• A * 

i 1 i , 1 1 " 





20 
15 

5 



0.10 




0.20 



FIG. 12. The amplitude of the Fourier peaks A,^, 
A2v, and the symmetry parameters a and b as a func- 
tion of anisotropy parameter |cr4| of the interfacial free en- 
ergy. The kinetic anisotropy has been varied proportion- 
ally {Pi = —204). Other relevant parameters are A = 0.55, 
v = m, ^ = 0.3 and ao = 0.16. 
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FIG. 13. Pattern formation under periodically alter- 
nating heating and cooling of amplitudes 60 = 8 (a) and 
64 (b). Other relevant parameters are A = 0.6, v = 60, 
^ = 0.5, (T4 = —0.12, /34 = 0.24, while the elapsed time t = 
0.22. 
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FIG. 14. The average velocity of the tip under periodic 
heating and cooling (open circles; parameters as for Fig. 
|l^ ), and for oscillatory pressure (open triangles; calculated 
with < p >— 0, ^ = 0.5, other parameters as for Fig. |2|(d)). 




FIG. 17. The effect of oscillatory pressure on the smec- 
tic-B dendrite growing into undercooled nematic phase (AT 
= 1.0 °C). (a) Without pressure oscillations; (b)-(f) with 
square wave pressure modulation of parameters pe = 2 bar, 
C = 0.2, u = 0.75, 0.91, 1.18, 1.64, and 1.96 Hz, respectively. 




FIG. 15. Pattern contours formed under periodic heat 
pulses of amplitudes increasing from left to right {bo — 5.7, 
6o = 28.6, and bo — 51.4). Other relevant parameters are 
A = 0.61, z/ = 50, 5 = 0.15, (74 = -0.16, /34 = 0.32, 
h = -14.3, To = 16 and At = 7 x 10"^ 
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FIG. 16. The average velocity Vo (open triangles) and 
the amplitude of velocity oscillation (full circles) as a 
function of the heating amplitude bo for simulations pre- 
sented in Fig. |l^. 
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FIG. 18. The effect of the v.o^...c.. 



tic-B dendrite growing into undercooled nematic phase 
(AT^ 1.0°C). (a) No oscillatory heating; (b) v = 0.2 
Hz, P = 6 X lO"'^ W/cm^ = 0.14; (c) i> = 0.46 
Hz, P = 6 X 10-' W/cm^, 5 = 0.16; (d) v = 0.81 Hz, 
P = 1.1 X 10"* W/cm^ ^ = 0.61. The two lines in panels 
(a) and (c) mark the positions, where the width w of the 
dendrite has been measured as a function of time (for the 
Fourier transforms see Fig. ^9|). Black lines in panel (b) 
denote the position of the tip at the centers of the heating 
pulses. 
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FIG. 19. Fourier spectrum W(f) of the half-width of the 
dendrite measured 66.5 /im behind the tip as shown in Fig. 
|l^(c). For comparison, the respective spectrum for free 
dendritic growth [Fig. b_q(a)] is also presented (see insert). 




FIG. 20. The effect of the amplitude of oscillatory pres- 
sure on the smectic-B dendrite growing into undercooled 
nematic phase (AT — 1.0 °C). (a) Without pressure os- 
cillations; and with square wave pressure modulation of 
parameters v = 1.64 Hz, ^ — 0.2, pe = 1 bar (b) and 2 bar 



(c), respectively. 
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FIG. 21. Side-branch formation close to the tip un- 
der "symmetric" (^ — 0.5) square-wave heating as dis- 
played by a smectic-B dendrite growing into the under- 
cooled nematic phase. The growth pattern (insert) and 
the Fourier-spectrum W{f) (main frame) are shown. Black 
lines in the insert denote the forcing frequency v. The white 
line marks the position (53.2pm behind the tip), where the 
width w of the dendrite has been measured as a function of 
time. Note the presence of the mode oscillating with twice 
the base frequency of forcing; a mode missing from the forc- 
ing spectrum. The relevant parameters are: i> = 0.35 Hz 
and P = 9 X lO"'' W/cm^ 
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